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ABSTRACT 

We study the inverse Compton scattering of solar photons by Galactic cosmic-ray electrons. We show 
that the 7-ray emission from this process is substantial with the maximum flux in the direction of the 
Sun; the angular distribution of the emission is broad. This previously-neglected foreground should 
be taken into account in studies of the diffuse Galactic and extragalactic 7-ray emission. Furthermore, 
observations by GLAST can be used to monitor the heliosphere and determine the electron spectrum 
as a function of position from distances as large as Saturn's orbit to close proximity of the Sun, thus 
enabling unique studies of solar modulation. This paves the way for the determination of other Galactic 
cosmic-ray species, primarily protons, near the solar surface which will lead to accurate predictions of 
7-rays from pp-interactions in the solar atmosphere. These albedo 7-rays will be observable by GLAST, 
allowing the study of deep atmospheric layers, magnetic field(s), and cosmic-ray cascade development. 
The latter is necessary to calculate the neutrino flux from pp-interactions at higher energies (>1 TeV). 
Although this flux is small, it is a "guaranteed flux" in contrast to other astrophysical sources of 
neutrinos, and may be detectable by km'^ neutrino telescopes of the near future, such as IceCube. 
Since the solar core is opaque for very high-energy neutrinos, directly studying the mass distribution 
of the solar core may thus be possible. 

Subject headings: elementary particles — Sun: general — Sun: interior — Sun: X-rays, gamma rays 
— cosmic rays — gamma-rays: theory 



1. INTRODUCTION 



Interactions of Galactic cosmic-ray (CR) nuclei with 
the solar atmosphere have been predicte d to be a source 
of ver y high energy (VHE) neutrinos (Moskalenko et al.l 
Il99lt iSeckel et al.l Il99ll: llngel man fc ThunmanI 119961) . 
and T-ravs (jSeckel et al.lll991|) . Thev are the decay prod- 
ucts of charged and neutral pions produced in interac- 
tions of CR nucleons with gas. The predictions for these 
albedo 7-rays give an integral flux F^{> 100 MeV) ~ 
(0.2 - 0.7) X lO"'^ cm-2 s^i, while analysis of EGRET 
7-ray telescope data has yielded only the upper limit 
F^{> 100 MeV) ^ 2.0 x lO"'^ cm'^ s'^ (jThompson et al.l 
1199 it) . At lower energies (<100 MeV) a contribution 
from CR electron bremsstrahlung in the solar atmo- 
sphere may exceed that from Tr'^-decay (see Fig. 4 in 
IStrong at al]|2004al ). 

Cosmic-ray electrons comprise ^1% of the total CR 
flux. However, they propagate over the heliospheric vol- 
ume, which is large compared to the solar atmosphere 
where the albedo 7-rays are produced. The Sun emits 
photons that are targets for inverse Compton (IC) scat- 
tering by CR electrons. As a result the heliosphere is a 
diffuse source of 7-rays with a broad angular distribution. 
In this paper, we evaluate the importance of IC scatter- 
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ing within the heliosphere, and discuss the consequences 
of its measurement by such instruments as the upcom- 
ing Gamma Ray Large Area Space Telescope (GLAST) 
mission. In the following, we use units h = c ~ = 1 . 

2. ANISOTROPIC INVERSE COMPTON SCATTERING IN 
THE HELIOSPHERE 

The distribution of CR electrons within the heliosphere 
is approximately isotropic. However, the distribution 
of solar photons is distinctly anisotropic, with the pho- 
tons propagating outward from the Sun. The expres- 
sion for the IC production spectrum for an arbitrary 
photon angular distr i bution is given by eq. (8) from 
iMoskalenko fc Strong! ()2000f ): 



di?(7e,ei 
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where is the classical electron radius, 7e is the elec- 
tron Lorentz-factor, ei and are the energies of the 
background photon in the laboratory system (LS) and 
electron rest system, correspondingly, £2 is the LS en- 
ergy of upscattered photon, £2 < 2"fef-'i/{'^ + 2e'j),e'j = 
ei7e(l + /?e cos C), C is the LS angle between the momenta 
of the electron and incoming photon (C = for head-on 
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Fig. 1. — Variables involved in the calculation. The large plot 
shows a 2-dimensional sketch; the dotted lines are a projection of 
the 3-dimensional angles. The inset defines the angles in 3 dimen- 
sions. 

collisions, see Fig. [J), £2 max = 4ei7g/(l + 4ei7e) is the 
maximal energy of the upscattered photons, and Qui^f) 
is the normalised angular distribution of target photons 
at a particular spatial point (/ dfl,^ Qi, = 1). 

The 7-ray flux as a function of 9 (see Fig. [Ij can be 
calculated by integrating eq. ([T]) along the line-of-sight 
taking into account the distribution of electrons and solar 
photons: 
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where the factor 1/4 comes from the angular distribu- 
tion of photons for the cas e of an emitting surface (see 
iMoskalenko fc Stronej|2000l eqs. [24]-[26]), Je(r, 7e) is the 
electron intensity, r is the radial distance from the Sun, 
Rq is the solar radius, and dnbb{ti,TQ) / dei is a black- 
body distribution at temperature Tq. 
Close to the Sun, Qv{^u) is given by 
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sinasin(/)COS ?/;, 

)y (r, 0) d cos (/> = 1, sin0™^'^ — 

Rq/t. For sufhciently large distances, this reduces to 
a delta function, Qu{d) — (5(cos C — cos a), where cos a = 
{d- acose)/{a'^ + - 2adcos6')i/^ a = 1 AU, and the 
variables are illustrated in Fig. (T] 

The IC 7-ray flux (eq. [2]) essentially falls as 1/r with 
heliocentric distance. Therefore, it is straightforward to 
estimate the region which contributes most to the IC 
emission. The inner boundary (in AU) is given by 



sin 9, 
1, 



9 < 90°, 
9 > 90°. 



Approximately 90% of the total 7-ray flux is produced 
between ri and r2 = lOri (Figure [2]). 

The electron spectrum as a function of position in the 
heliosphere is obtained by adopting the force-field ap- 
proximation and deriving the radial dependence of the 
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Fig. 2. — The range of heliospheric distances probed by the IC 
7-rays versus angular distance from the Sun. 90% of the predicted 
7-ray flux is produced in the region between the solid and dashed 
lines. Dotted lines show the solar size. 



modulation potential using a toy model. This should 
be sufficient, as we seek only to illustrate the effect, 
and are interested in electron energies above ~5 GeV 
where the heliospheric modulation is moderate. Assum- 
ing a separability of the heliospheric diffusion coefficient 
K = 0Ki( r,t)K2(l,t), the rnodula tion potential can be 
obtained (jGleeson fc Axfordl[l96l : 



$(r, t) 



3Ki(a:, t) ' 



(6) 



where t is the time variable, rb is the heliospheric bound- 
a ry, and V is the solar win d speed. 

iFuiii fc McDonaidI ()2005[ ) have analysed the radial de- 
pendence of the mean free path. A, during the solar min- 
ima of Cycles 20/22 and Cycle 21 using the data from 
the IMP 8, Voyagers 1/2, and Pioneer 10 spacecraft. For 
Cycles 20/22 they obtain a parameterisation A oc r^-^ in 
the outer heliosphere and A oc r^'^ in the inner helio- 
sphere, where the break is at 10-15 AU. For Cycle 21, 
an r^'^ dependence fits well for both the inner and outer 
heliosphere. Using these parameterisations, ki oc A, and 
assuming V = const yields (Cycles 20/22) 
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where $o is the modulation potential at 1 AU, ro = 
10 AU, rb — 100 AU, and we have neglected the time 
dependence. For Cycle 21 we have 
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3. CALCULATIONS 



The local interstellar CR electron spectrum is ap- 
proximated by dJe(rfc,ee)/dee = 2 x 105(ee + 400)-3-22 
cm^^ s^^ sr^^ MeV~^, where Cg is the electron energy 
in MeV. This expression matches the interstellar local 
electron spectrum above ~500 MeV as calculated by 
GALPROP in plain diffusion and reacceleration models 
(jPtuskin et al.l 12006.) . The s olar photosph eric tempera- 
ture is taken as T = 5770 K (jLivingstonl l 19991. 

Figure [3] shows the integral intensity above 100 MeV 
and 1 GeV versus the angular distance from the Sun 
for no modulation, and for modulated electron spectra 
according to eqs. ([7]) and ^ with two modulation lev- 
els $0 = 500, 1000 MV; these correspond approximately 
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Fig. 3. — Integral intensity above 100 MeV and 1 GeV. Black 
line: no modulation; red lines: potential $1; blue lines: potential 
<l>2. Dashed lines: $0 = 500 MV; dotted lines: $0 = 1000 MV. 

to the solar minimum and maximum conditions, respec- 
tively. The dependency on $0 is strong and thus the 
emission will vary over the solar cycle while the depen- 
dency on the modulation potential parameterisation is 
much weaker. 

The 68% containment radius of the EGRET point 
spread function is ~ 6° at 100 MeV. For 6* < 6°, we 
calculate an integral flux F-,{> 100 MeV) 1.6 x lO"'^ 
cm~^ s~^ (modula ted). The sum of the IC and albedo 
emission Meckel ctHlllQQl over this region is consis- 
tent with the EGRET upper Umit of 2 x lO"'^ cm^^ g-i 
(jThompson et al.lll997l ). 

The differential intensities for selected angles 6 are 
shown in Fig. ID The parameterisations $1.2, eqs. ([7]) 
and ([8|), yield very similar differential 7-ray intensities; 
we thus show the results for <i>i only. Also shown is 
the extragala c tic 7-ra y background (EGRB) obtained by 
IStrong at al.l (|2004af) . which includes true diffuse emis- 
sion and an unresolved source component. For 9 < 5°, 
the solar IC emission is more intense than the EGRB. 
The latter is the major component of t he diffuse 7- 
ray em ission outside of the Galactic plane (jStrong at al.l 
l2004b( l. Table [1] gives the integral intensities above 10 
MeV, 100 MeV, and 1 GeV averaged over the whole sky. 
The values for 500 MV and 1000 MV are obtained by av- 
eraging over the fluxes for $1 and $2- Above 100 MeV, 
the sky-a veraged intensity is ~ 5-10% of the EGRB ob- 
tained bv IStrong at all ()2004aD . 

Note that unresolved blazars are thought to be the ma- 
jor contributors to the EGR B (jStecker &: Salamoiilll996t 
iKneiske fc MannheimI [20 05) ; thous ands of them will be 
resolved by GLAST. Therefore, the IStrong at all (|2004af) 
EGRB flux is an upper limit for the true diffuse emission. 

4. PERSPECTIVES FOR GLAST 

The Large Area Telescope (LAT) high-energy 7-ray 
detector, in preparation for launch by NASA in late 
2007, will have unprecedented angular resolution, ef- 
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Fig. 4. — Differential intensities for selected 8. Line-sets (top to 
bottom): 0.3°, 1°, 5°, 10°, 45°, and 180°. Solid line: no modu- 
lation; dasfied line: <I>o = 500 MV; dotted l ine: 'I'o = 1000 MV . 
Data points: diffuse extragalactic 7-ray ffux llStrong at al.ll2"004al ). 



fecti ve area, and fiel d of view for high-energy 7-rays 
(Mc Enerv et al.ll2004D . The LAT will scan the sky con- 
tinuously and provide essentially complete sky coverage 
every 2 orbits (approximately 3 hours). As a conse- 
quence, the light curve for any detectable source will 
be well sampled. Also, the scanning motion will permit 
frequent evaluation of the stability of the performance 
of the LAT with respect to "standard candle" sources. 
Based on the expected sensitivity of the LAT'^, a source 
with flux ~10~^ cm^^ s^^ and the hardness of the solar 
IC emission will be detectable on a daily basis when the 
Sun is not close to the Galactic plane, where the diffuse 
emission is brightest. Sensitive variability studies of the 
bright core of the IC emission surrounding the Sun should 
be possible on weekly time scales. With exposure accu- 
mulated over several months, the Sun should be resolved 
as an extended source and potentially its IC, Tr'^-decay, 
and bremsstrahlung components separated spatially. 

The extended IC emission is not isotropic, but accu- 
mulated over year timescales will be uniform in ecliptic 
longitude and brightest at low ecliptic latitudes. The 
ecliptic plane crosses the Galactic plane near the Galac- 
tic Centre, and the solar IC component may be impor- 
tant for investigating the nature of t he reported "halo" 
of 7-rays about the Galactic Centre (jPixon et al.| [l998). 

The solar IC contribution to the overall celestial diffuse 
emission can be modelled directly as we have done here. 
Measurement of the spectrum of solar IC emission in 
the near Sun direction fixes the emission at all angles. 
This is especially true at high energies where the solar 
modulation in the outer heliosphere is negligible. Also, 
for any particular direction of the sky the Galactic and 
extragalactic components can be derived from the data 
by measuring variations in diffuse intensity with solar 
elongation angle. 

5. DISCUSSION AND CONCLUSION 

In this paper we have studied the IC scattering of solar 
photons by CR electrons'^. We have shown that the emis- 
sion is significant and broadly distributed with maximum 

http://glast.stanford.edu 
^ Wh en this work had already b een completed we learned about 
work bv lOrlando fc Strong! (200^ on the same subject. 
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Table 1. All-sky average integral intensity 



E 


So = 


500 MV 


1000 MV 


>10 MeV 


5.9 


3.7 


2.5 


>100 MeV 


0.7 


0.6 


0.5 


>1 GeV 


0.05 


0.05 


0.05 



Note. — Units 10~* cm"^ sr~ 



brightness in the direction of the Sun. The whole sky is 
shining in 7-rays contributing to a foreground that would 
otherwise be ascribed to the Galactic and extragalactic 
diffuse emission. 

The IC emission from CR electrons depends on their 
spectrum in the heliosphere and varies with the modu- 
lation level. Observations in different directions can be 
used to determine the electron spectrum at different he- 
liocentric distances. A sensitive 7-ray telescope on orbit 
could monitor the heliosphere, providing information on 
its dynamics. Such observations also could be used to 
study the electron spectrum in close proximity to the 
Sun, unreachable for direct measurements by spacecraft. 
The assumed isotropy of the electron distribution every- 
where in the heliosphere is an approximation. Closer to 
the Sun, the magnetic field and non-isotropic solar wind 
speed affect the CR electron spectrum and angular distri- 
bution, which in turn will produce asymmetries in the IC 
emission. Observations of such asymmetries may provide 
us with information about the magnetic field and solar 
wind speed at different heliolatitudes including far away 
from the ecliptic. 

What are the immediate implications? Sin ce solar 
modulation theory is well develo ped fe.g.. Florins ki et all 
l2003HFerreira fc iPotgiete^l2004D . accurate measurement 
of the electron spectrum near the solar surface will open 
the way to derive the spectra of other Galactic CR 
species, primarily protons, near the Sun. The CR pro- 
ton spectrum is the input to calculations of 7-rays from 
pp-interactions in the solar atmosphere. These predic- 
tions can be further tested using GLAST observations of 
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pionic 7-rays from the Sun. In turn, this would provide 
information about the density profile of the solar atmo- 
sphere, magnetic fi eld(s), and CR cascade development 
(jSeckel et al.l 119911 ). The higher the energy of 7-rays, 
the higher the energy of the ambient particles, and thus 
the depth of the layers tested is increased. In conjunction 
with other solar monitors this can bring understanding of 
the deep atmospheric layers. Sun spots, magnetic storms, 
and other solar activity. 

Furthermore, understanding the solar atmosphere 
is necessary to calculate the neutrino flux from pp- 
interactions at higher energies (>1 TeV), as the CR 
shower development depends on the density distribu- 
tion and the und erl ying magnetic field structur e (see 
ISeckel et allll99 1' and'Ingel man fc Thuninaiilll996l for de- 
tails). Although small, the VHE neutrino flux from the 
Sun is a "guaranteed flux" in contrast to other astro- 
physical sources of neutrinos. It may be detectable by 
km'^ neutri no telescopes of the near future, such a s Ice- 
Cube (e.g. JHettlage et al.ll2000HAhrens et"alll2004[ ). on a 
timescale of a year or several years. Therefore, only long 
term periodicities associated with the solar activity could 
be detecta ble. Since the sola. r core is opaque for VHE 
neutrinos (ISeckel et al.l 119911: [Moskalenko fc Karakulal 
II993I : llngelman fc ThunmanI 11996! ). observations of the 
neutrino flux may provide us with information about the 
solar mass distribution. 

To summarise, we have shown that the observation of 
the IC emission from CR electrons in the heliosphere 
and its distribution on the sky will open a new chapter 
in astrophysics. This makes a sensitive 7-ray telescope 
a useful tool to monitor the heliosphere and heliospheric 
propagation of CR, and to study CR interactions in the 
solar atmosphere. 
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ERRATUM: "INVERSE COMPTON SCATTERING ON SOLAR PHOTONS, HELIOSPHERIC MODULATION, 

AND NEUTRINO ASTROPHYSICS" (ApJ, 652, L65 [2006]) 

Igor V. Moskalenko, Troy A. Porter, Seth W. Digel 



We noticed an error in the description of the distribution of solar photons at an arbitrary distance from the Sun, 
equations (3) and (4). The correct expression is 

-1 

(1) 

fl--fj <cos(/)<l, (2) 

i.e. Qi,(r, (p) is independent of (j) within the soUd angle covered by the Sun. Applying the correct angular distribution 
does not give results that are noticeably different from those obtained with the delta-function (pure radial) photon 
distribution. Indeed, it should be the case since in the energy range under consideration 7e ^ 1 and the ambient 
photon angular distribution can be approximated by the delta-function. 

We also discovered a numerical error in the code which affects the results below ^1 GeV, especially in case of small 
9. Figures [3] and m show the corrected integral and differential intensities. Table 1 shows the corrected all-sky average 
integral intensities. The 68% containment radius of the EGRET point spread function is ~6° at 100 MeV. For 9 < 6°, 
the corrected integral flux is F^{> 100 MeV) ^ (2.0 — 4.3) x 10"^ cm~^ s~^, where the given range corresponds to 
different modulation levels ($o = 1000 - 500 MV). 
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Fig. 3. — Integral intensity above 100 MeV and 1 GeV. Black 
line, no modulation; red lines, potential <I>i; blue lines, potential 
4>2. Dashed lines, <I>o = 500 MV; dotted lines, <I>o = 1000 MV. 
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Fig. 4. — Differential intensities for selected 8. Line-sets (top to 
bottom): 0.3°, 1°, 5°, 10°, 45°, and 180°. Solid line: no modu- 
lation; dashed line: *o = 500 MV; dotted line: $o = 1000 MV. 
Data points: diffuse extragalactic 7-ray flux (A. W. Strong, I. V. 
Moskalenko, & O. Reimer [ApJ, 613,956 (2004)]). 



Table 1 

All-sky average integral intensity 



E 


-I'D = 


500 MV 


1000 MV 


>10 MeV 


7.1 


6.5 


6.0 


>100 MeV 


1.3 


1.2 


1.1 


>1 GeV 


0.14 


0.13 


0.11 



Note. — Units 10 ^ cm ^ s ^ sr 



